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We review a representation of Hubbard-like models that is based on auxiliary pseudospin variables. 
These pseudospins refer to the local charge modulo two in the original model and display a local 
Z2 gauge freedom. We discuss the associated mean-field theory in a variety of different contexts 
which are related to the problem of the interaction-driven metal-insulator transition at half-filling 
including Fermi surface deformation and spectral features beyond the local approximation. Notably, 
on the mean-field level, the Hubbard bands are derived from the excitations of an Ising model in a 
transverse field and the quantum critical point of this model is identified with the Brinkman-Rice 
criticality of the almost localized Fermi liquid state. Non-local correlations are included using a 
cluster mean-field approximation and the Schwinger boson theory for the auxiliary quantum Ising 
model. 
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I. INTRODUCTION 

Strong correlation physics is a central ingredient 
for diverse solid-state systems including the high-T c 
cuprates, 1 other transition metal oxides 2 or unconven- 
tional superconductors 3 as well as fractional quantum 
Hall states. 4 Related physics is currently also discussed in 
the context of ultracold atoms in optical lattices. 5 Strong 
interactions can give rise to a variety of unusual quan- 
tum phases including ordered phases in spin, charge, and 
orbital degrees of freedom, as well as miscellaneous ex- 
otic liquid phases. 6 The complexity of these many-body 
systems lies in the fundamental dichotomy between real 
space (localization) and momentum space (derealiza- 
tion) in combination with the restrictions on the Hilbcrt 
space which are enforced by the strong interaction among 
the particles. Only a few rigorous results are available for 
comparable "simple" models such as the Hubbard or the 
t-J Hamiltonian in more than one 7,8 or less than infinite 
dimensions. 9,10 Hence, there is a considerable amount of 
ongoing work attempting to reach the physically most 
relevant limit of two or three dimensions. 11 Thereby, the 
classical problem of the Mott metal-insulator transition 12 
plays a key role and its nature in various systems is still 
actively debated. Recent examples include the orbital- 
selective Mott transition in multiorbital systems such 
as Ca2- a; Sr :c Ru04, 13 the problem of momentum space 
differentiation related to the pseudogap phenomena in 
cuprates 14 or the paramagnetic metal-insulator (spin liq- 
uid) transition in frustrated geometries such as in the 
organic compound k-(BEDT-TTF) 2 Cu 2 (CN) 3 . 15 

Widely used theoretical approaches to tackle these 
problems are slave-particle methods because they pro- 
vide powerful tools to deal with the restrictions im- 
posed on the Hilbert space due to strong correlations. 
They have been pioneered in the context of quantum 
magnets, 16,17 magnetic impurities in metals 18-20 and 
doped Mott insulators. 21 ' 22 The basic idea is to repre- 
sent local degrees of freedom with the help of auxiliary 



degrees of freedom in an enlarged Hilbcrt space. How- 
ever, in order to have a faithful representation, these aux- 
iliary degrees of freedom obey certain constraints and are 
not independent (although in some cases they survive as 
"real" particles) 1 but are "enslaved" - hence, the name. 
In practice, the starting point is usually a mean-field 
state in the enlarged Hilbert space obeying a set of self- 
consistency equations. This allows for a semi-analytical 
and non-perturbative treatment of correlation effects. 

In this article, we make contact with a particular 
class of slave-particle representations for Hubbard-like 
models. They share a simple physical picture of the 
interaction-driven paramagnetic metal-insulator transi- 
tion which dates back to early works of Gutzwillcr 23 and 
Brinkman and Rice 2 and others 25 who introduced the 
notion of an "almost localized Fermi liquid" to char- 
acterize the metallic state close to the Mott transition. 
In particular, the transition to the localized (insulating) 
state is signaled by a diverging effective mass. In infi- 
nite dimensions, this picture can be put on a more firm 
ground 10 since correlations are strictly local but in two 
or three dimensions non-local correlations can lead to dif- 
ferent conclusions. Our approach is closely related to the 
four-boson method introduced by Kotliar and Rucken- 
stein (KR), 26 and its extensions. 27-34 In these approaches 
it is, however, not straightforward to include fluctuations 
of the mean fields and the high-energy (incoherent) part 
of the single-particle spectrum. Most studies are there- 
fore restricted to the low-energy (coherent) part of the 
one-particle spectrum (but note also Refs. [35-39]). An 
elegant formulation which can overcome some of these 
shortcomings has been given by Florens and Georges 40,41 
in terms of a slave-rotor representation. In this formula- 
tion, the phase fluctuations of the rotors give rise to the 
incoherent spectral features in the single-particle spec- 
trum. Moreover, it also yields a closer connection to the 
superfluid to insulator transition in the Bosc-Hubbard 
model. 42 We also note the slave-spin representation of 
dc'Mcdici and coworkers 43 which has been introduced to 



study the orbital-selective Mott transition in a two-band 
Hubbard model and follows a similar spirit. 

We contribute to these different approaches by review- 
ing an alternative slave-spin formulation which has re- 
cently been applied to the study of dynamically generated 
double occupancy in cold atomic Fermi systems. 44 The 
advantage of our formulation is that it reduces the com- 
plexity of the representation to a minimum. The gauge 
freedom is only Z2 and the auxiliary quantum model, 
which describes the high-energy degrees of freedom in the 
mean-field approximation, is given by the transverse-field 
Ising model. Furthermore, in the simplest treatment, we 
exactly recover the result of the Gutzwiller approxima- 
tion applied to the metal-insulator transition. 

The outline is as follows: we first introduce the gen- 
eral formulation of the problem using slave-pseudospin 
variables in Sec. II. We then discuss the mean-field ap- 
proximation which consists of decoupling pseudospin and 
fermion degrees of freedom in Sec. III. Furthermore, 
in Sees. IV and V we investigate the consequences of 
two subsequent approximations made for the pseudospin 
problem: (i) the single-site mean-field approximation and 
(ii) the use of Schwinger bosons to treat the fluctua- 
tions around the (renormalized) classical ground state. 
In Sec. VI we discuss the importance to include an aver- 
aged local constraint. 



II. SLAVE-SPIN FORMULATION 

We begin this section by introducing the general frame- 
work of the slave-spin formulation we want to utilize in 
the study of the Hubbard model, Eq. (1). We define a 
representation of physical operators in an enlarged local 
Hilbert space and we analyze the additional local sym- 
metry which is introduced by this procedure. The subse- 
quent discussion of the non-interacting case allows us to 
set the stage for the mean-field study in the remainder 
of the paper. We conclude this section by commenting 
on Elitzur's theorem (stating the impossibility to break 
the aformentioned local symmetry) and the restrictions 
it poses on the interpretation of the mean-field results 
discussed later. 

For concreteness, we shall consider the single band 
Hubbard model written in the form 



U 



(i) 



The hopping amplitude between sites i and j is denoted 
by tij and U is the onsitc repulsion. cj£? destroys (cre- 
ates) an electron at site i with spin a and hi — c L c io-- 
Throughout the paper we work at half filling, thus assum- 
ing (hi) = 1. 

We now introduce a representation of the local physical 
states which is based on auxiliary pseudospin variables, 
see Fig. 1. Thus, let us introduce a pseudospin I with 
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FIG. 1: (color online) Pictorial illustration of the pseudospin 
representation discussed in this article. 
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encoding doubly occupied and empty sites (|+)) and 
singly occupied sites (|— )). Consequently, the eigenvalue 
of I z refers to the presence (—1/2) or absence (+1/2) of a 
local magnetic moment. In addition, auxiliary Fermi cre- 
ation and annihilation operators f„ are introduced to 
preserve the canonical anti-commutation relations (see 
App. A for a connection to earlier work). The physical 
creation (annihilation) operator of the original model is 
then represented as 



(3) 



The physical states in the enlarged local Hilbert-space 
are 

|e) = |+)|0), \p a ) = \-)\a), |d) = |+)|2), (4) 

where |0) is the vacuum of the /-fermions, 

|<r)=/t|0) and |2) = /t/]|0). 

The states (4) are pictorially shown in Fig. 1. In the lat- 
tice system the above definitions are generalized for each 
lattice site i. The physical subspace is selected by the 
requirement that the local quasiparticle charge modulo 
two can be represented by If as follows: 



(m - i) 2 = 0. 



(5) 



Associated with the constraint (5) , let us define an oper- 
ator 



\-(m- i) 2 



1 



If [l-2(n 4 -l) 2 ] 



Qi has eigenvalues and 1 and the eigensectors define 
the local physical subspace and its orthogonal com- 
plement respectively. Thus, the local Hilbert space 
Hi is decomposed according to 
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The projection onto the physical subspace can then also 
be achieved by imposing for each lattice site 

Qi = 0, V*. (6) 

As a result, we can write the Hubbard interaction in the 
physical subspace solely using pseudospin operators If 
and the original Hamiltonian is represented as 

h> = -4 £ /, ; ///;.c./; ; , + \ £ (z? + 1) • (7) 

As long as the constraint (6) is fulfilled, the Hamiltonian 
(7) is equivalent to the original model (f). In the fol- 
lowing we often drop the constant term UN s /2 (N s — 
number of sites) in Eq. (7). Using Qi, the projector onto 
the physical subspace takes the form 

V = l[(l-Q l ), (8) 

i 

with V 2 = V and VH'V is equivalent to the original 
model (1). 

A. Gauge structure 

In the following section, we will briefly discuss the 
gauge freedom introduced by the representation (4). Al- 
though the gauge structure is not important for the 
mean-field analysis outlined in the remainder of this ar- 
ticle, it allows us to gain a deeper understanding of the 
slave-spin construction. 

The gauge group is formed by all local transforma- 
tions which leave the local physical annihilation and 
creation operators cfj — 21 f f^J (and therefore also 
H') invariant. We will see that the most obvious Z 2 - 
transformations 

c la = 2I?f ia = 2{-I?){-f ia ), 

cl = 2Iffl = 2(-7f )(-/*,), (9) 

are in fact not the most general ones but are part of a 
larger U(l) group. Nevertheless, as shown below, the Z 2 - 
gauge transformations (9) are the only ones respecting 
the mean-field decomposition and in this sense, the mean- 
field ansatz breaks the U(I) gauge symmetry down to Z 2 . 

In order to explicitly derive all possible gauge trans- 
formations we make use of the fact that H' does not mix 
the physical subspace with its complement. The gener- 
alized local charge Qi is therefore a conserved quantity 
and commutes with H': 

[Qi,H'} = 0. 

A conserved charge Qi arises from a continuous sym- 
metry Ui(4>i) = exjp(icf>i Qi) with < 4>i < 2n. In 



other words, Qi generates the local gauge transforma- 
tions Ui(cf>i) and we conclude that there is a local U(l) 
freedom, 

^(^0^(0= cfj, (10) 

with c$ = Under the action of U^fa) the 

pseudo-fermions transform as 

Ui(<f>i)£>Uifc)1 =e i ^£\l-Q i ) + e- i ^f^Q i , (11) 

and the pseudospins as 

UifaWUitftf = e^Ifll-QJ+e-^IfQi, 
UifaWUifa? = If. (12) 

The transformation relations (11) and (12) allow to ex- 
plicitly check the relation (10). 

At this point it is important to recall that the present 
slave-spin scheme is a projective construction: a physi- 
cal state |\&) is obtained from a general state |<&) of the 
enlarged Hilbert space after projection, |W) = V\<&). A 
different state |$') which is obtained from |$) by a gauge 
transformation will give rise to the same physical state 
l^). (For a related statement in the projective construc- 
tion of spin liquid phases see e.g. Refs. 45 and 46.) This 
is an intuitive way to understand the origin of the gauge 
freedom. In practice, we start from a mean-field (prod- 
uct) state in pseudospin and pseudo-fermion degrees of 
freedom: |$mf) = |$/)|$/), see Sec. III. Applying the 
transformations (11) and (12) to the mean-field Hamil- 
tonians (19) and (21) we find that in general only gauge 
transformations with fa = or n respect the mean-field 
product form. This means that different mean-field states 
which give rise to the same physical state are related by 
gauge transformations with = or 7r. In this sense, 
the mean-field ansatz breaks the U(l) freedom down to 
the smaller Z 2 freedom. Therefore, we anticipate that 
the relevant gauge group to study fluctuations around 
the mean-field state is Z 2 . It consists of the following 
two elements: the identity operator, id^ = Ui(<fii = 0), 
and the operator 

Ui = Ui{4n = tt) = (-1) Q * =\-2Q t . (13) 
The action of U{ on f^J and If is given by 



ujfut = -If. 
Obviously, Ui promotes the gauge transformations (9). 

B. Non-interacting model U = 

We now turn to the discussion of the non-interacting 
model: 

//; i^ /,,/;• /;/;;./•„. 
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The eigenstates of Hq are (at least in principle) exactly 
known. This fact allows us to study the action of the pro- 
jector V and the role of the local gauge freedom in a more 
rigorous manner. Moreover, the non-interacting limit is 
helpful for the interpretation of the mean-field approxi- 
mation introduced and explored later in this article. In 
the following we use arguments which are similar to those 
given in a discussion of an exactly solvable spin-model on 
the square lattice in Ref. 47. 

Since [If, Jig] = we choose the eigenstates |<E>) 
of H (in the enlarged Hilbert space) as being prod- 
uct states in the (Ising-) pseudospin and fermion de- 
grees of freedom. Explicitly, |$) = | {cx}} | <?i({o})) where 
If \{a}} = a.i\{a}), on — ±1/2. The pseudo fermion com- 
ponent |</>({a})) is a Slater-Determinant obtained from 
the non-interacting model with effective hopping param- 
eters AaitijCtj. As an example, let us choose oti = +1/2 
for all sites i. The ground state in this sector is given 
by |$ ) = |{l/2})|0„({l/2})> where |^({l/2})) is the 
Fermi sea of the pseudo fermions. Clearly, |$o) bas the 
same energy Eq as the physical ground state: 

However, the state |$o) is not an eigenstate of the pro- 
jector V and, therefore, is not the physical ground state. 
Instead, the (unnormalized) physical ground state is ob- 
tained by projecting |4>o) onto the physical subspace, 

At this point it is important to note that the local Z2 
gauge freedom (9), [H , uf\ — 0, implies that each energy 
sector of H is macroscopically degenerate (~ 2^*). In 
fact, applying any combination of local gauge transfor- 
mations Ui to |$o) yields an eigenstate of Hq with energy 
Eq but which is (in general) orthogonal to |$o)- This 
property allows to obtain the physical ground state by a 
particular superposition of eigenstates in the -Eo-sector. 
To see this, let us take a closer look at the projector. 
Using Eq. (13) the projector Eq. (8) is written in the 
form 

i 

Explicitly, it takes the form 

V= fl + ^ Ui + u n u l2 + ■ ■ ■ +\\uA /2 N \ 

\ i ii<l2 i / 

(15) 

The action of Ui on an eigenstate |$) of Hq can be un- 
derstood by writing in the form 

Ui = -2J?[l-2(ni-l) 2 ]. (16) 

Therefore, u$ changes the sign of Qfj in the pseudospin 
component |{ck})- To understand the action of Ui on the 
pseudo fermion part, we expand |</>({a})) m the site di- 
agonal occupation number basis. The form (16) of Uj, 
implies that the components with an empty or doubly 



occupied site i are multiplied by —1 while those with 
a singly occupied site i are not changed. The resulting 
wave-function is then just the corresponding eigenstate 
of a non-interacting model in which f^J is replaced by 

— fit ■ Thus, Ui indeed acts as a local gauge transforma- 
tion: 

u l (\{a})ma}))) = -\{a'})ma'})), (17) 

where {a 1 } = {. . . , ctj-i, — ot^ oti+i, . . . }. Hence, ow- 
ing to the form (15) of V , the physical ground state 
|*o) = V\Qq) is the equal amplitude superposition of 
all the degenerate states which are obtained from |$o) 
by applying all possible local gauge transformations. 
The physical subspace is therefore the gauge invariant 
subspace. 

Note that, although |$q) IS not the physical ground 
state, it is a representative state of the ground-state en- 
ergy sector and, as long as gauge invariant operators 
O = UiOui are considered, the expectation values are 
equal, ($ |O|$ > = (*o|0|* )- 49 Because physical ob- 
servables are gauge invariant we can calculate all physi- 
cal properties of the non-interacting model by restricting 
to the sector on — 1/2 for all i (or to any other fixed 
configuration {a}). 

C. Relevance of Elitzur's theorem 

In this paragraph, we briefly comment on Elitzur's 
theorem 50 which states the impossibility to sponta- 
neously break a local symmetry. This implies that ther- 
mal averages (. . . )th in the physical as well as in the 
enlarged Hilbert space of operators which are not gauge 
invariant have to be zero. In the present case, it implies, 
for example, that (If )th = 0. 51 

In the non-interacting limit discussed above this result 
is plausible. Let us first consider thermal averages in the 
enlarged Hilbert space: because each energy sector of H' 
is spanned by a macroscopic number of states involving 
all possible pseudospin configurations {a} the expecta- 
tion value of If averages to zero. Likewise, restricting 
the thermal average to the physical subspace, the fact 
that physical states are equal amplitude superpositions 
of all the degenerate states in the enlarged Hilbert space 
results in a mutual cancellation of positive and negative 
contributions to (If)th- In particular, for the physical 
ground state we find (* |-f?|*o) = 0. 

However, the expectation value of a single state in the 
enlarged Hilbert space can have a non-vanishing expec- 
tation value {^f\If |$) 7^ 0. Consider for example the 
exact eigenstate |$o) of Hq defined in the previous sec- 
tion. Although (&o\If |$o) = 1/2 we can correctly ob- 
tain all the physical observables in the ground state from 
|$o) because it is a representative state of the physical 
ground-state sector. 

We would like to argue that the mean-field approxima- 
tion introduced in the next section should be interpreted 
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in the same sense. Namely, in the mean-field approxi- 
mation, we seek for a product state which approximates 
one particular ground state of H' in the enlarged Hilbert 
space. In fact, for small U, the mean-field state is contin- 
uously connected to |$o)- While this procedure does not 
yield a systematic approximation we assume that also 
for larger U the mean-field state is sufficiently close to 
a representative state of the ground-state energy sector. 
Moreover, the mean-field approximation is truly varia- 
tional: it gives an upper limit for the true ground state 
energy of the original model. Although a mean-field state 
will in general have (If) ^ 0, we can use it to approxi- 
matively calculate physical ground state properties. 

We think that a similar interpretation of other slave- 
particle mean-field theories is appropriate. This would 
be consistent with observations made by studying slave- 
boson theories in the "radial gauge" where all the ex- 
act physical properties of a toy model have been ob- 
tained from the saddle-point solution of a functional 
integral. 52,53 



III. MEAN-FIELD THEORY 

Let us now discuss the mean-field theory which re- 
sults from the representation introduced in the previ- 
ous section. To this end, we assume product states 
\9) = \^fj)\^ff) in pseudospin and fermion degrees of 
freedom. These states live in the enlarged Hilbert space. 
As discussed in the previous section, we assume that \\fr) 
yields a sufficiently good approximation to a state in the 
(macroscopically degenerate) ground-state sector of H' 
in the enlarged Hilbert space. We postpone a discussion 
of the importance to include the constraint (5) on average 
to Sec. VI. Note, however, that the relation 

W> + 5 = ^W-((E^-i) a ), (18) 

holds, which leads to the identification of the physi- 
cal fraction of doubly occupied sites, (clfC.^cl^c^) , with 
(If) /2+1/ 4. This is analogous to earlier mean-field treat- 
ments, see Eq. (A3). The rational behind the mean-field 
decoupling is the fact that we can approximately distin- 
guish two energy scales. Indeed, as shown below, the 
characteristic scale of the pseudospins is ~ max(U, U c ) 
whereas that of the pseudo-fermions is ~ t [U c is given 
in Eq. (27)]. Although this observation together with 
the insights gained from the non-interacting limit justi- 
fies to some extent the mean-field decoupling, it should be 
considered as a first step on which a more sophisticated 
analysis can be based. Nevertheless, on the present level 
of approximations we can make close contact to earlier 
results. 



A. Mean-field Hamiltonians 

As a consequence of the mean-field decoupling we ob- 
tain two effective Hamiltonians: the fermion problem as- 
sumes the form of a non-interacting tight-binding Hamil- 
tonian, 

H f = (*/|fT'|*i) = - J2 mUiflJiv, (19) 

with the hopping amplitude Uj renormalized by a factor 

9ij = i(lfl*)j. (20) 

On the other hand, the pseudospin problem reduces to 
the quantum Ising model 

H Z = (9 f \H'\* f ) = - £ xnUjlfl? + \ £ If, (21) 

(id) « 

with the transverse field U/2 and the Ising exchange cou- 
pling UjXtj where 

Xij = + c.c). (22) 

a 

The sum over in Eq. (21) means summation over 

all bonds where Xij 7^ 0. 

The quantum Ising model (21) is a prime example 
of a system displaying a quantum critical point at a 
critical ratio of the transverse field to the Ising cou- 
pling. It separates a magnetically ordered from a quan- 
tum paramagnet. 54 Below we show that this quantum 
critical point can be identified with the Brinkman-Rice 
criticality of an almost localized Fermi liquid. 24,25 

Equations (20) and (22) are the two coupled self- 
consistency equations to be solved in the mean-field ap- 
proximation. We note that there is always the trivial 
solution gij = Xij = of these equations which is the 
physical solution in the atomic limit £y = 0. Non-trivial 
solutions invoke an (approximate) solution of the quan- 
tum Ising model. 

The symmetry properties of the original model are 
conserved in the mean-field approximations (19) und 
(21). For example, particle number conservation leads 
to a global U(l) symmetry for the slave fermion sector 
[Eq. (19)] in the usual manner. Within our approach, 
the pseudospin sector explicitly breaks the pseudospin- 
rotation symmetry [Eq. (21)]. This is not related to the 
particle number, however, because the pseudospins mea- 
sure charge only modulo 2. Therefore, the appearance of 
terms cx Ifl? (instead of I^Ij) do not bias our system 
toward a U(l) symmetry broken phase in the mean-field 
approximat ion . 

B. Single-particle Green function 

Eventually, we are interested in physical (gauge- 
invariant) quantities such as the single-particle Green 



G 



function, 



G„(v i ,r j ;t) = -i{Tc j(r (t)cl(0)), 



where T denotes time ordering. In the mean-field theory, 
G a is obtained as 

G ff (r i ,r i ;t) Ri 4B ij (t)G{(r i) r J -;t) ) (23) 

where we have introduced the auxiliary quantities 

By(t) = (Tif(t)if(O)), 

G^,r i; t) = -i{Tf j(T (t)fl(0)). 

In momentum and energy space, the relation (23) 
translates into a convolution of £?(q, ui) and G£(q, w). 
It is noteworthy to mention that the canonical anti- 
commutation relations of the physical annihilation and 
creation operators are preserved on average, 

({Cfar, cL}> = 4<Jf if )<{/ iff , 4,}) = %<W, (24) 



where (...) denotes the average over mean- field eigen- 
states. As a consequence, the single-particle spectral 
density is correctly normalized, as long as the spin-1/2 
identity (if) = 1/4 is respected. 



IV. MEAN-FIELD APPROXIMATIONS TO THE 
ISING MODEL 

A straight forward way to study the quantum Ising 
model (21) is mean-field approximations. In the simplest 
case we consider a single pseudospin coupled to a self- 
consistent effective field. This is the local approximation 
discussed in Sec. IV A. In order to improve over the local 
approximation we consider in Sees. IV B and IV C the 
pseudospin problem on a finite cluster; see Refs. 55 and 
56 for related works. This allows us to discuss important 
aspects of inter-site correlations which are absent in the 
local approximation. In the following we work in the 
zero temperature limit and we restrict our analysis to 
translation-invariant and paramagnetic states. 



A. Local approximation 

We start with the single-site "cluster". Noteworthy, 
on this level of the approximation, the mean-field self- 
consistency (20) and (22) leads to the Brinkman-Rice 
transition 24 obtained in the paramagnetic Gutzwiller ap- 
proximation of the Hubbard model. 23 To see this, let us 
introduce the mean magnetization and 



H 



MF 



= hY.it, h=^ Jut + 



jx\2 



(25) 



where the pseudospin has been rotated due to the action 
of the mean field U c (Iq), 



Ii = e iaI he- iaI i , tana 



(26) 



In Eqs. (25) and (26) we have used the dimensionless 
interaction parameter u = U/U c where the energy scale 
U c is associated with the pseudospin Ising coupling in the 
single-site solution 

U C = Y *tf *j = ~ 16 / dee Pa (e) > 0. (27) 

J-D 



Here, /0<j(e) is the non-interacting density of states per 
spin with half bandwidth D and ep is the effective Fermi 
energy of the pseudo fermions. 

Self-consistency of (Iq) yields the pseudospin magne- 
tization and from Eq. (20) we obtain the hopping renor- 
malization factor as follows: 



4(i, 



x\2 




(28) 



In particular, on this level of the approximation, the ef- 
fective mass m* /to = I j g diverges at the critical interac- 
tion strength U c , indicating the transition to the localized 
state. The double occupancy follows from the relation 
(18): 



(IS) 



0. 



u < 1; 
u > 1. 



(29) 



The angle a of the pseudospin rotation (26) can be writ- 
ten in terms of d as 



cos - = \/\ - 2d 2 . 



(30) 



Following Eq. (30), the rotation of the pseudospin cor- 
responds to adjusting the average fraction d 2 of doubly- 
occupied sites. 

Note that the single-site solution (28) and (29) for 
u > 1 reproduces the result of the atomic limit tij = 0. 
Due to the fact that inter-site correlations have been ne- 
glected the low-energy physics of the Mott insulator is 
completely absent. This shortcoming can be addressed 
by going beyond the local approximation. 

Our mean-field state has a finite expectation value 
(if) 7^ for u < 1. However, according to Elitzure's 
theorem, the physical ground state requires (If) — be- 
cause if is not a gauge invariant operator. As we have 
argued in Sec. II C we consider our mean-field state as an 
approximation to a representative eigenstate in the en- 
larged Hilbert space. In this light, we can still expect that 
reasonable approximations for physical (gauge invariant) 
observables are obtained. 



B. Cluster approximations 

To improve over the local approximation we have stud- 
ied clusters with two or more sites. We briefly review the 
generic differences. A particular example is discussed in 
the next subsection. Most important is the fact that in 
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finite spatial dimensions the value of (if) 2 and (If If) 
with i ^ j are different. In other words, inter-site cor- 
relations are introduced. Consequently, in accordance 
with previous studies, 41 ' 56 there is a distinction between 
the quasiparticle weight Z and the effective mass renor- 
malization m/m* of the quasiparticles. While in gen- 
eral the factor g = m/m* stays finite across the Mott 
transition, the quasiparticle weight Z still vanishes for 
U — >• U c . The quasiparticle weight and the effective mass 
in the mean-field theory follow from the form (23) of the 
single-particle Green's function. For a nearest-neighbor 
hopping model they are given by 



0.15 



Z = 4(Jf ) 2 



and 



= g = mm 



for nearest-neighbor pairs i,j. The distinction between 
g and Z is now apparent in finite dimensions. More gen- 
eral, the electronic self-energy obtains a k dependence 
which, in the metallic phase at particle-hole symmetry, 
is of the form 



E(w,k) = (1 + Z~ 1 )u + 



9_ 

Z 



l)£k 



for k near the Fermi surface defined by £k = and for 
small w's. This issue will again be discussed in Sec. V 
when we consider the Schwinger boson mean-field the- 
ory in order to access the low lying excitations of the 
quantum Ising model. 

In addition, the critical interaction strength for the 
Mott transition is rcnormalized compared to the value of 
the local approximation (27). We find, however, that the 
exact value depends on the choice of the cluster. 



C. Deformation of Fermi surface 

As an application of the cluster mean-field scheme, let 
us now study the question how inter-site correlations 
can change the shape of the Fermi surface. Such an 
interaction-driven deformation is easily described within 
the present scheme when considering a generic disper- 
sion with further-neighbor hopping amplitudes. The de- 
formation then results from a different renormalization 
of nonequivalent hopping amplitudes. As an example, 
we study here a model on a two-dimensional square lat- 
tice with nearest-neighbor hopping amplitude t and next- 
nearest-neighbor hopping amplitude if . The renormal- 
ized quasiparticle dispersion is then given by 

£k = —2g t t (cos k x + cos k y ) — Agt't' cos k x cos k y — ep. 



Here, the renormalization factors are < 



4(Jf If) with 



■ 9 t C) -\'t- 3 , 
i and j (next-)nearest neighbors. The parameter is de- 
termined to satisfy the Landau-Luttinger sum rule for the 
Fermi surface defined by £k = 0. According to Eq. (22), 
the auxiliary quantum Ising model acquires the nearest 
neighbor exchange 

Xt = 4i J^(cos k x + cos fc z )n kcr , 

k,<r 



3.25 




0.05 



3.15. 



FIG. 2: The nearest-neighbor exchange \t an d the next- 
nearest neighbor exchange yfi! °f the auxiliary pseudospin 
model as function of U/t for t' = —0.3t. The inset shows the 
considered cluster. 



and a next-nearest neighbor coupling 

x't' = 8t' cos k x cos fc 9 nkff. 

We have considered a 2 x 2 cluster as shown in the inset 
of Fig. 2 and have solved the self-consistency Eqs. (20) 
and (22) for different values of U/t with a fixed ratio 
t'/t = —0.3. Figure 2 shows \ an d x' as function of U/t. 
Note the difference in scale for \ an d x' ■ Figure 3(a) 
shows the quasiparticle weight Z along with g t and gt>. 
As mentioned in the previous subsection, Z vanishes in 
the Mott insulator while gt and g# stay finite. On the 
metallic side, g t > < gt which means that the Fermi surface 
is deformed toward the fully nested surface with t' = 0, 
as shown in Fig. 3(b). A similar behavior has been found 
in calculations which take antiferromagnetic fluctuations 
into account. 57 On the insulating side, g# gt- This can 
be expected since gv/gt ~ J'/J ~ (t'/t) 2 in the Mott 
insulator where J (J 1 ) is the super-exchange for (next-) 
nearest neighbors. For U > U c , the low-lying excitations 
are not given by Landau quasiparticles since Z = 0. The 
surface for U — 20t shown in panel (b) can be viewed as 
the Fermi surface of a "hidden Fermi liquid" 58 or as the 
spinon Fermi surface of a gapless spin liquid. 1 However, 
since this Fermi surface is close to perfect nesting we 
expect that (residual) interactions open a gap. Before we 
study the effect of such residual interactions in Sev. VI 
in more detail, we present another way of going beyond 
the local approximation in the next section. 



V. SCHWINGER BOSON THEORY 

An alternative approach to go beyond the local approx- 
imation is the use of the Schwinger boson theory for the 
quantum Ising model and we follow an essentially similar 
line of thoughts as in Ref. 59. The bosonic representation 
establishes also a connection to earlier slave-boson repre- 
sentations, see App. A. Moreover, it allows us to take the 




FIG. 3: (color online) (a) The quasiparticle weight Z and the 
nearest-neighbor and next-nearest-neighbor hopping renor- 
malization factors g t and g t i , respectively, as function of U/t. 
(b) The Fermi surface for different values of the interaction 
strength. 

effect of quantum fluctuations into account and it yields 
the low-lying excitations. 

We introduce two sets of Bose creation and annihila- 
tion operators xf^ and to represent the pseudospin 
algebra 

In order to obtain a faithful representation of the spin 
algebra, these operators have to obey 

In the Schwinger boson formulation, the pseudospin ro- 
tation given in Eq. (26) translates to an unitary transfor- 
mation of the Bose creation and annihilation operators 




The a and b bosons can be interpreted as the Schwinger 
bosons of the rotated pseudospin I,. Equally, using the 
relation (30), we can specify the transformation (31) by 
the value of d and in the following we denote the canon- 
ically transformed Ising model by Hg(d). 

A. Classical ground state 

The result of the local approximation is re-obtained by 
assuming a product form of the wave function 

|O>=n«ll0), (32) 

i 

and optimizing the energy 

E(d) = (0\H B (d)\0), 

with respect to the parameter d. We recover again the 
Gutzwiller result (28). The state (32) is the classical 



ground state of the transverse Ising model. We state 
here that this is the point where in the KR functional 
integral representation one stops if only the saddle point 
(without fluctuations around it) is considered. 

B. Role of fluctuations in three dimensions 

The formalism developed so far offers a framework to 
study the excitation spectrum of the transverse Ising 
model. We show that these excitations lead to the in- 
coherent one-particle excitations of the Hubbard model 
(upper and lower Hubbard bands). Furthermore, quan- 
tum fluctuations renormalize the classical ground state 
obtained in the mean-field approximation. In general, 
the role of fluctuations crucially depends on the dimen- 
sionality of the system. Here we restrict our analysis to 
the three dimensional (cubic) lattice where fluctuations 
around the classical ground state are small for most pa- 
rameters. Nevertheless, the finite dimensionality reintro- 
duced in our analysis causes interesting changes in the 
nature of the Mott transition as compared to the infinite 
dimensional result. Our results are in agreement with 
previous results based on the slave-rotor formalism. 41 

1. Effective Hamiltonian for fluctuations 

Formally, the effective Hamiltonian is derived by ex- 
panding Hs{d) up to second order in the b bosons. The 
parameter d of the unitary transformation is then deter- 
mined by the requirement that quadratic mixing terms in 
a and b of Hb id) vanish. This yields again the condition 
(29). The next step is to let the a bosons condense and 
to replace the a operators by unity. In momentum space 
the effective Hamiltonian then reads 

H-^^r^h <33) 

for u < 1. For u > 1 we find 

HeS = ^y / Bl(^ + U M W (34) 

We have introduced the operators given by 

4 = (&L*-k) and Sk= (fet k k )' 

as well as 

1 3 

7k = --^cos(fcj). 
6 i=i 

The effective Hamiltonian for fluctuations, Eq. (33) and 
(34), is diagonalized by the following Bogoliubov trans- 
formation 

f &k \ /coshtf k sinhtfk \ / /3 k \ 

V b -k J \ Sinll,9 k COsht? k ) { /?l k ) ■ 
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FIG. 4: (a) The band of pseudospin-wave excitations obtained 
in the spin- wave analysis of the transverse Ising model in three 
dimensions as function of u — U /U c . (b) The inverse effective 
mass m/m* and the quasiparticle weight Z as function of 
U /U c - Note that at the Mott transition u — 1, Z vanishes in 
contrast to m/m*. 



The mixing angle is given by 



•die — -atanh 
2 



min(l,M 2 )7 k 



min(l, it 2 )7k + 2 max(l, u) 



The /3-operators describe the low- lying eigenmodes of the 
transverse field Ising model and correspond to (gapped) 
pseudospin-wave excitations. The spectrum of these aux- 
iliary excitations and their relation to physical properties 
is discussed in the next paragraph. 



2. Pseudo-spin-wave mode and Mott-Hubbard gap 

From the diagonalization of the effective Hamiltonian 
we find the following pseudospin-wave spectrum 



^ = Uc x f y /l + u 2 7k , for u < 1, 
2 \ \/u 2 + U7k, for u > 1, 



(35) 



with an excitation gap A = hu>(0). The quantum critical- 
ity at u — 1 is reflected in the softening of the mode (35). 
For u > 1, the jump A/i in the chemical potential from 
hole to particle doping amounts to twice the excitation 
gap, 



An = 2A = [Ml 



1 



(36) 



The above pseudospin-mode corresponds to the gapped 
charge excitation of the Mott insulator and Eq. (36) co- 
incides with the expression found for the Mott-Hubbard 
gap in the Kotliar-Ruckenstein formulation. 35,37 The 
band of the pseudospin mode (35) is shown in Fig. 4(a) 
for different values of the interaction strength u. 



3. Renormalized ground state 



find 

Eg 
N s 



u,. 



(1 



BZ 



47T 3 



(l - v/l+^Tk) , 



BZ 



d 3 k 
4tt 



| (l - V 7 ! + 7k/«) , 



where the integrals over the Brillouin zone (BZ) represent 
quantum corrections to the result of the Gutzwiller ap- 
proximation. Note that Eq < for any finite u, meaning 
that the trivial mean-field solution with g = \ = is al- 
ways higher in energy for any finite u. This is in contrast 
to the single-site mean-field approximation where the so- 
lution for u > 1 coincides with the atomic limit. The 
present scheme based on the Schwinger boson representa- 
tion of the pseudospin is analogous to taking into account 
Gaussian fluctuations around the mean-field transition 
without renormalizing the actual transition point. 



Effective mass and guasiparticle weight in three 
dimensions 



As in the cluster mean-field treatment above and in 
previous studies invoking slave rotors, 41 we obtain that 
fluctuations reintroduce inter-site correlations beyond 
the mean-field value. In particular, there is a distinction 
between the U dependence of the quasiparticle weight 
and of the effective mass renormalization of the quasipar- 
ticles, see Fig. 4(b). While the hopping renormalization 
factor g — m/m* stays finite across the Mott transition, 
the quasiparticle weight Z still vanishes for u — > 1. In 
the insulating phase, the effective mass is obtained from 
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irr' 



de p a (s) 



e/D 



1 + -e/D 



u, 



2UJ dep < ,(e)(e/D)* = ^ 7 , 

where D — 6t is half of the band width and the second 
line holds in the limit a> 1. Thus, similar to the cluster 
mean-field calculation, the pseudospin-wave analysis cap- 
tures the energy scale of the superexchange J = 4t 2 /U 
for large U and m/m* — J%/(8i) is finite due to vir- 
tual hopping processes which keep the number of doubly 
occupied sites small but finite. It is interesting to note 
that this value is the same as in the uniform resonating- 
valence-bond (w-RVB) phase of the Heisenberg model 
Hg = </X)(i j) S< "Sj which is the strong-coupling limit of 
the low-energy sector of the Hubbard model at half fill- 
ing. The u-RVB phase is obtained by a uniform spinon 
mean-field ansatz for a single real Hubbard-Stratonovich 
field x [given by Eq. (22)] i 1 



o 5^ ifiafj 



h.C. 



Quantum fluctuations lead to a renormalization of the 
ground-state energy Eq (u < 1) and Eq (u > 1). We 



which yields the same effective mass (w/m*) u _ RVB — 
Jx/(8t). The transition from a metal to a gapless spin 
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liquid in frustrated geometries has recently been dis- 
cussed by several authors. 60-66 



5. Spectral one-particle density 

We now turn to the discussion of how the auxil- 
iary pseudspin excitations affect the physical excitation 
spectrum. To this end we calculate the spectral one- 
particle density A a {uS). We start by writing A a (uj) in 
the Lehmann representation 



+ (0|4» 6{uj + uj n0 ) 



where |n) denotes an eigenstate of the full Hamiltonian 
with energy E n and ui nm = E n — E rn . In the slave-spin 



method used here, the true eigenstates are approximated 
by the mean-field eigenstates obtained in the pseudospin- 
wave analysis. We then have to calculate matrix elements 



of the form (0|ig7, 



The details of the calculation are 



lengthy but straightforward, 
results. 



Here we present the main 



The spectral weight contains a coherent quasiparti- 
cle sector A™ h (w) as well as an incoherent contribution 
A™ c (w). We find that the coherent contribution is given 
by 



= - Pa (oj/g). 
9 

Note that A a (0) = A c f h (0) cx Z/g gradually vanishes 
when approaching the Mott insulator, in contrast to the 
infinite dimensional result. In the metallic phase u < 1, 
the incoherent contribution is dominated by 



4?» 



AD 



c JQ 



Pa 

de p a {e) x { p<y 
0, 




A_ < ui < A + + gD; 
-A + - gD < uj < A_; 
else; 



(37) 



where A± — Hf\/l±u 2 denote the edges of the excitation spectrum. On the other hand, in the insulating phase 
u > 1, the coherent part vanishes ^4™ h (u;) = and we find 



A a (ui) = A^ c {uj) 



4D 



de Pa {e) x <^ Pa 
0, 



D I i(u-g£) 

» 1 u% 

D f 4(u+ge) 2 

u { U? 



A_ < w < A + + gD; 
-A + - gD < u < A_ 
else; 



(38) 



where A± = ¥ y/l ± lju. The spectral density is 
shown in Fig. 5 for different values of the interaction 
strength. The gapped mode (35) found in the trans- 
verse Ising model leads to the incoherent weight around 
± max(J7 c , U)/2 in the spectral density. 37 In the metal- 
lic phase, we find the characteristic three peak struc- 
ture with preformed Hubbard bands centered at hui ~ 
±Z7 c /2 and a coherent Gutzwiller band at hui m 0. The 
Gutzwiller band disappears at u = 1 and the Hubbard 
bands touch at hui — 0. From the expression (38) we find 
that in the large U limit the Hubbard bands assume a 
constant width of Whb = U c /2 s» 8t and are separated 
by U. This is consistent with the large U expansion of 
the jump in the chemical potential (36) and the general 
expectation that both an added hole or an added double 
occupancy are mobile with an energy — Whb/2: 



U-Whb. 



The reduction in the bandwidth Whb of the upper (or 
lower) Hubbard band as compared to the noninteracting 
bandwidth W = 2D = 12t is in qualitative agreement 
with the retraceable-path approximation for a single hole 
doped into an infinite-C/ Mott insulator. 67 



6. One-particle sum rule and fluctuation regime 

To estimate the validity of the mean field plus spin- 
wave calculation, we compare the fluctuations with the 
magnitude of the order parameter and define u R through 
the condition 



1 



1 + 6 



(CM) 



= (In* 



For the cubic lattice considered here we obtain u R w 0.9, 
hence, for a«110.1 the fluctuation induced corrections 
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to the mean-field result are important and the validity 
of the above analysis is limited. In particular, the ap- 
proximative nature of our treatment of the pseudospin 
problem violates (If) 2 — 1/4 and the spectral weight 
fails to be properly normalized. As shown in Fig. 6, the 
one-particle sum rule 



duAa-fjjj) = 1, 



is fulfilled within 10%-12%. The failure predominantly 
manifests itself in the fluctuation regime, however. Vi- 
olations of sum rules are a known shortcoming of the 
Gaussian approximation. We expect that some of these 
inconsistencies can be cured by taking into account ap- 
propriate mode-mode couplings. 



^ 1.1 

b 

3 

^ 1 



0.9. 



Fluctuation 
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FIG. 6: The integrated spectral weight of the one-particle 
spectral density in the pseudospin- wave analysis. The shaded 
region denotes the fluctuation regime where the magnitude of 
the fluctuations is of the same order than the order parameter. 



VI. ROLE OF THE CONSTRAINT 

So far, we have discussed various aspects of the mean- 
field approximation where the eigenfunctions |^) = 
of the slave-spin Hamiltonian (7) in the en- 
larged Hilbert space have been approximated by prod- 
uct states in pseudospin and pseudo fermion degrees of 
freedom and we have not paid attention on the relation 
(5). There is an obvious shortcoming in this treatment: 
the expectation value (7f) + 1/2 is in general not equal 
to ((n, — l) 2 ) where n, = fLfia 1S tne local pseudo- 
fermion density. But since both operators are gauge in- 
variant and the expectation values are equal for the phys- 
ical state by construction we conclude that this relation 
should in fact hold for all eigenstates |$ m ) of if': 



1) 2 |$, 



(39) 



Therefore, we expect to obtain a better approximation 
of an eigenstate of H' by enforcing the relation (39) for 
mean-field states. This can be achieved by the usual La- 
grange multiplier method. In this way, we can access 
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FIG. 5: (color online) The contour plot shows the one-particle 
spectral density A a (uS) with the coherent Gutzwiller and the 
(preformed) Hubbard bands as function of energy u/t and in- 
teraction strength u. On the right-hand side, A a (uj) is shown 
for fixed values of it in the metallic and the insulating phase. 



a larger class of correlations in the mean-field approxi- 
mation and a prominent example missed so far is mag- 
netism. In the following we discuss two possibilities to 
satisfy (39) for mean-field states. The first approach re- 
lies on the presence of a finite staggered magnetization v s 
and is therefore expected to be of relevance for magneti- 
cally unfrustrated models. The second approach involves 
an iteration of the slave-spin approximation in the para- 
magnetic sector and its natural applications are therefore 
models where magnetic order is fully frustrated. In this 
brief overview, we do not analyze weakly frustrated sys- 
tems such as the t — t' model discussed in Sec. IV C. For 
these models the situation is more complex since aspects 
such as the Fermi surface deformation are likely to play 
a role as well. 

We follow a standard procedure to enforce the relation 
(39) for mean-field states by adding a term 



A 



■l5> 



It + 2 - ("* - !) 2 



(40) 



to the slave-spin Hamiltonian (7). One then seeks for a 
stationary solution with respect to the Lagrange multi- 
plier A^: 



d(H' + A) 
dX,. 



(41) 



A. Neel state in unfrustrated systems 

On a bipartite lattice a Neel-ordered state is usually 
a leading instability and in the following we explore 
the possibility that the term (40) triggers this instabil- 
ity. Furthermore, introducing an order parameter for 
the staggered magnetization v s allows to satisfy the con- 
straint (5) on average, see also Eq. (39). Hence, we de- 
couple (40) in the antiferromagnetic channel by writing 
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(at half filling) 



-y ("it - ^4-) + 2 4 A 



where z/, 



'4, 



and we assume V; 



~> (42) 
if i belongs 



to the A sublattice and v t = —v s for i on the B sublattice. 
Furthermore, we seek for a translational invariant state 
with Aj = A. In addition to the slave-spin self-consistency 
(20) and (22) there are two more self-consistency equa- 
tions in order to determine the parameters g, x, A and 
v.: 
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X 

2 
1 
A 



mil) 



-</?), 

r sp=Q 



2 1 de- — 
- zt y/Ag'e 



2 C 2 



(43) 
(44) 



Equation (43) is just the average constraint using the de- 
coupling (42) and (44) is the gap equation for the pseudo- 
fermion problem. 

We have solved these equations for different values of 
U/t on a square lattice with nearest-neighbor hopping t. 
The result for v s is summarized in Fig. 7 for the local 
slave-spin approximation (Sec. IV A) and with fluctua- 
tions (Sec. V B) . Also shown is the result obtained in the 
standard Hartree-Fock (HF) mean-field theory applied 
to the original model (1). We find a first-order tran- 
sition in the staggered magnetization both in the local 
approximation and with the inclusion of fluctuations in 
the pseudospin-wave sector. Moreover, for a finite range 
of U /t, the self-consistency equations have two distinct 
solutions for v a which both locally minimize the energy. 
One solution corresponds to a "low-zV phase of delocal- 
ized character ("Slater") the other to a "high-zAj" phase 
of localized character ( "Mott-Heisenberg" ) . (However, 
the naive decoupling (42) should be taken with care, in 
particular for larger values of U/t.) 

As seen in Fig. 7, the present approach coincides with 
the pure HF of the original model only in lowest order 
in U/t, but clearly differs in higher orders. These pro- 
nounced differences are attributed to a "positive feedback 
effect" between the staggered magnetization v s and the 
hopping renormalization factor g which, formally, enters 
by the averaged constraint (43). In the \<w-v s phase, the 
system optimizes the kinetic energy by sustaining a weak 
staggered magnetization but a small effective mass (large 
g). v s is therefore smaller than in HF. On the other hand, 
the high-iAj solution results from the optimization of the 
potential energy: the physical mechanism behind it is the 
fact that an enhanced effective mass (localization) favors 
the formation of local magnetic moments (singly occu- 
pied sites) which allows to minimize the potential energy 
cost. v s is therefore larger compared with the HF result. 

Including fluctuations in the slave-spin problem re- 
duces the regime where the high-t/ s solution is stable. It 
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FIG. 7: (color online) The staggered magnetization u s at 
T = for the square lattice as function of the interaction 
strength U /t calculated using the local slave-spin approxima- 
tion of Sec. IVA(black, dashed-dotted) and with fluctuation 
corrections as discussed in Sec. V B (red, full). For comparison 
we also show the result obtained in the Hartree-Fock mean- 
field theory (blue, dashed). In the slave-spin calculation there 
is a finite range of U/t where there exists two solutions for v a 
which corresponds to a (local) minima of the energy. 



is likely that the quantum fluctuations in v s , which have 
not been considered here, will reduce the stability regime 
even more and it is an open problem if the first-order 
transition survives once we go beyond the mean-field de- 
coupling (42). (We note that in dynamical mean-field 
theory on the infinite dimensional hypercubic lattice a 
crossover is observed.) 68 Furthermore, such fluctuations 
lead to an overall reduction in the staggered moment 
also in the large U/t limit. From spin- wave theory 69 
and Monte Carlo simulations 70 of the spin-1/2 Heisen- 
berg model on the square lattice it is known that the 
sublattice magnetization is only around 60% of its clas- 
sical value. 



B. Iterated slave-spin scheme for full frustration 

We now discuss an alternative way to satisfy the con- 
straint on average which does not rely on the presence 
of a magnetic order parameter. This scheme is therefore 
more likely to play a role in fully frustrated systems where 
magnetic order is suppressed. Also in this approach we 
find a first-order transition between an itinerant and a 
localized phase as function of U/t. 

The term (40) introduces an onsite interaction of 
strength A (we assume A.; = A) for the / fermions in 
the mean-field treatment and Hf takes the form of a 
Hubbard model with renormalized parameters: effective 
hopping gitij and onsite interaction A. Let us again use 
the slave-spin mean-field theory to treat this renormal- 
ized model: following Sec. II we introduce an additional 
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set of pseudofermions h\ and pseudopins K; to repre- 
sent the / fermions and the eigenstates of Hf are ap- 
proximated as product states \^ff) = |^E'^)|^'/ l ). In order 
to satisfy the original constraint between the / fermions 
and the I spins on average, we can relax the constraint 
between the h fermions and the K spins. The reason is 
that the relation (18) for the (h,K.) pair gives 

(7 

But from (41) it follows that (K-) = (If) and the above 
relation is just the average constraint for the (/, I) pair. 
Equation (41) also implies that the ratio of the transverse 
field to the exchange coupling is the same for the I and 
the K model and that gj = gx = g. As a result, wc 
find that A = {7/2. In the local approximation for the 
transverse field Ising model, a non-trivial solution exists 
only for v < 1 where v = 4u/(3-\/3). The kinetic energy 
is renormalized by g 2 where for v < 1 

l + 2cos(0/3) 
.9= 3 , 

with = arccos(l — 2v 2 ) and g — otherwise. We thus 
find a first order transition at a reduced critical interac- 
tion strength 4t/ c /(3v3) between a paramagnetic metal 
and an insulator. 

We note here that a paramagnetic first-order transition 
at T = has recently been reported 66 for the single-band 
Hubbard model on the (highly-frustrated) triangular lat- 
tice. These calculations revealed a first-order transition 
from a paramagnetic metal to a paramagnetic insulator. 
Moreover, a second first-order transition at larger val- 
ues of U/t between the paramagnetic metal and the 120° 
Neel-ordered state was found. 



VII. CONCLUSIONS 

We have reviewed a slave-particle formulation for 
strongly interacting Fermi systems which can be con- 
sidered as a minimal formulation of previous represen- 
tations. Our approach is based on auxiliary pseudospin 
variables displaying a local Z 2 gauge freedom. They re- 
fer to the local charge modulo two. The simplicity of the 
representation allows to exactly solve the non-interacting 
model and to gain insights into the artificial local sym- 
metry introduced by enlarging the Hilbert space. For 
the interacting model, the main focus so far has been the 
investigation of the corresponding mean-field approxima- 
tion. We have proposed an interpretation of the mean- 
field decoupling in terms of approximating a "represen- 
tative eigenstate" in the enlarged Hilbert space rather 
than the physical eigenstate. In this light, our mean- 
field results do not contradict Elitzur's theorem and we 
have argued that reasonable results for gauge invariant 
observables can be obtained. 



The mean-field theory has been applied to various as- 
pects related to the interaction-driven Mott transition in 
the single-band Hubbard model. The role of inter-site 
correlations present in finite dimensions has been dis- 
cussed in a cluster approximation and we have found that 
inter-site correlations deform the Fermi surface toward 
the fully nested one. By including quantum fluctuations 
in the auxiliary Ising model, we have analytically calcu- 
lated the single-particle density of state and have elabo- 
rated the connection between the (gapped) pseudospin- 
wave excitations and the (preformed) Hubbard bands of 
the original model. 

We have also commented on the importance to enforce 
an averaged constraint which has to be satisfied by all 
eigenstates in the enlarged Hilbert space. In this way it 
is also possible to include magnetic correlations. For the 
nearest-neighbor hopping model on the square lattice we 
have found a first order transition between a phase with 
a low and high value of the staggered magnetization. 

There are many more open problems left for future 
studies and we hope to stimulate some further investiga- 
tions. For example, an extension of the present approach 
away from half-filling is highly desirable; in particular in 
view of potential applications of the current formalism 
to inhomogeneous systems (such as artificially structured 
correlated materials 33,34 or cold atomic gases in optical 
traps). Furthermore, we encourage an investigation of 
possible "flux patterns" (and its excitations) associated 
with the local Z 2 gauge freedom on the mean-field level. 
In order to go beyond the mean-field description it is an 
interesting problem to study the Z 2 lattice gauge the- 
ory coupled to the "matter field" of the pseudo fermions 
as obtained from the present approach. Moreover, such 
a study could also help to make contact with a previ- 
ously considered Z 2 gauge theory for strongly correlated 
electrons 22 and could add much to our understanding 
of the present formulation and its connection to various 
other slave-particle representations. 
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Appendix A: Relation to the four boson and the 
two-spin formulation 

In this appendix we relate our approach to the KR four 
boson formulation for a one-band model. 26 An equivalent 
representation is also obtained by using two slave spins 
introduced by deMedici and co-workers in Ref. 43. We 
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show that our representation is obtained by restricting 
the action of the two aforementioned slave spins to the 
particle-hole and spin-cr symmetric subspace. 

In the following, it is sufficient to focus on the local 
Hilbcrt space. Kotliar and Ruckenstein introduced four 
Bose creation and annihilation operators e^, p& and 
as well as two Fermi operators f„' to represent the 
local occupation number states as 



|e> = et|fi>, \a)=pifi\n), 



where |f2) is a fake vacuum state. The constraints to 
project out unphysical states are given by 



et e + 5jptp ff + t *td = 1, (Al) 

dU + pl Pa = flf a . (A2) 

The above constraints also imply that 

dU = /j/t/jA, (A3) 

which is analog to Eq. (5). The physical annihilation and 
creation operators can be represented as 



Ca = Zafa, 4 = 4/J, (A4) 

where the simplest form for the z-operators is 

z a = p\d + e)p a , 4 = d)p s +p\e. 

(In order to reproduce the correct non-interacting limit in 
the mean- field theory, a modified form was introduced.) 26 
The constraint (Al) can be resolved by introducing slave- 
spin variables 43 according to 

S£=zl, S~= z °i S z a + \=d)d + p\p IJ . 



The constraint (A2) then reads 



Su + 2 — ftrfa- 



(A5) 



We note that the slave-spin representation of Ref. 43 
differs from the one we have used in this paper since 
S* +1/2 refers to the presence (5| = +1/2) or absence 
(S* = —1/2) of a spin-cr electron. In contrast, in our 
method, only one slave-spin I is introduced and the I z - 
valuc refers to the total local charge modulo two. In the 
two slave-spin method, the physical annihilation operator 
was represented as 

c a = [S~ +a(n)S+]U, (A6) 

which in the physical subspace is equivalent to Eq. (A4) . 
The parameter 



a(n) 



- 1 



V n ( 2 ~ «) 

was adjusted such that the correct U = limit at particle 
density n is recovered. 56 In particular, for n = 1, the 
physical annihilation operator was written in the form 



(A7) 



Since (2S^) 2 = 1 it follows that for n = 1 the canonical 
anti-commutation relation is preserved in the mean-field 
decoupling, similar to Eq. (24). However, for n ^ 1, this 
property is lost. It seems that it is generally difficult to 
find a representation which preserves the anticommuta- 
tion relations and reproduces the correct noninteracting 
limit away from half-filling. 

The connection to our reduced slave-spin representa- 
tion at n = 1 is achieved by restricting the slave-spin 
operators to the two-dimensional particle-hole and spin- 
cr symmetric subspace spanned by 



10,1) 



1 

71' 



1,1 



i,i)). 



Here, we have introduced the triplet states \m z , 1) (m z — 
0, ±1) of the total spin S — Sf + S4.. The states |±) can 
be viewed as the eigenstates of a pseudospin operator I z , 
see Eq. (2). This establishes the connection between the 
different representations. 
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